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I.  Introduction. 

The  problem  to  be  considered  here,  that  of  a plane 
pressure  wave  impinging  on  a thin  spherical  shell,  was  sug- 
gested by  G.  F.  Garrier  in  consequence  of  work  previously 
done  by  him  on  a related  problem  [ 1]  • (cf.  also  [6].) 

An  attempt  had  been  made  to  determine  the  response 
to  an  incident  acoustic  wave  of  a thin  elastic  shell,  in  particu- 
lar a cylindrical  shell,  taking  into  account  both  the  incident  and 
diffracted  waves.  The  form  of  the  functions  dealt  with  in  the 
analysis  made  it  difficult  to  obtain  accurate  explicit  results. 

If  the  obstacle  is  taken  to  be  spherical  in  shape,  we  still  have 
a fairly  practical  though  highly  simplified  model  of  an  actual 
physical  structure;  moreover,  the  problem  becomes  mathematically 
simpler,  admitting  of  exact  solutions  for  the  deformation  and 
accompanying  strains  in  the  elastic  body.  It  was  therefore 
decided  to  treat  the  case  of  the  sphere  in  detail, 

1,  The  results  presented  in  this  paper  were  obtained  in  the  courss 
of  research  sponsored  by  the  Office  of  Naval  Research  under 
Contract  N7onr  - 35810  with  Brown  University, 

2.  Research  Associate,  Brown  University, 
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We  will,  as  for  the  cylinder,  deal  with  the  linearized 
theory  of  wave  propagation  in  a compressible  fluid,  and  with 
small  deflections  of  the  shell, 

II,  Forced  Vibrations  of  A Thin  Spherical  Shell , 

We  consider  a closed  shell  of  thickness  h with  h<<R 
where  R is  the  radius  of  the  middle  surface.  The  motion  of  any 
closed  oval*' shell,  in  particular  a spherical  shell  is,  by  a 
theorem  of  Jellett  [2],  primarily  extensional.  Therefore  the 
general  membrane  theory  of  shells  is  applicable.  If  we  locate 
the  origin  of  our  coordinate  system  at  the  center  of  the  shell, 
and  choose  as  the  z-axis  the  direction  of  propagation  of  the 
incoming  wave,  then  we  have  the  additional  simplification  of 
symmetrical  loading  (c.f.  Fig,  1).  The  equations  of  dynamic 
equilibrium  for  an  element  of  shell  may  therefore  be  written  [3] 

2 

JL(Nq  R sin  6)  - N R cos9  - ph  L-J  R2  sin0  = 0 (2,1) 

99  ™ at2 

2 

Nq  + N = (s  - ph  — ~)R  = 0.  (2,2) 

y 9t2 

Here  Nq,  and  Nq  are  the  normal  forces/unit  length  acting  on  the 
sides  of  the  element,  s is  the  applied  force  (radial  in  direc- 
tion), p is  the  shell  density,  and  v and  w are  the  tangential  and 
radial  components  of  the  displacement. 

To  eliminate  N^  and  Nq  from  e qua t ions  (2,1)  and  (2,2) 
we  make  use  of  Hooke’s  Law 

1,  Principle  radii  of  curvature  finite  and  of  the  same  sign. 
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e6  9 

e<pcp 


<N0  -vN9 
(N^  - v N0) 


(2.3) 

(2.k) 


and  the  expressions  for  the  strains  in  spherical  coordinates 

ee9  = SCf+^  (2*« 

e = - [v  cot  9 + w],  (2,6) 

TT  R 

Combining  (2,3),  (2,*+),  (2.5)  and  (2.6),  we  obtain  and  Nq 
in  terms  of  the  displacements  w and  v 


V Ne  = S;[|00t9  + 2I+if]  (2’7) 

Ne.Nf,.^qcote  -Ig]  (2.8) 

so  that  the  equations  of  equilibrium  become 


-S!_[v(-v+  m + aV  8W+  ax  ctn9  - V otn2  9 ] 
1-v2  99  a02  ae  99 

- ph  R2  = 0 

et2 

“~~[v  ctn  9 + 2w+  ■—  ] - [s  - ph  ]r2  = 0. 

1-v  99  r at2 


(2.9) 

(2.10) 


An  equation  in  v only  can  be  easily  obtained  as  follows. 
Consider  the  operator 

«o  ' 2 + a'^— ■ • <2.U) 

Mq(w)  is  found  from  (2.10)  to  be  equal  to 

ctn  0 - & + 1=1  R2s. 


-v 


(2,12) 
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If  Mq  is  then  applied  to  (2.9),  and  (2.12)  substituted  wherever 
Mq(w)  appears,  a relation  in  v only  results: 

M0C2.9)  = iz  + a^vf.ia-vip2^  eS  + 

0 E at2  e2  at^  VV 

+ v(l~v)p  K2  9^v  9^v  (l-v)pR2  a^v 

E 9t2  ” 2 90 2 E at2902 

+ 2v  ctn20  + L1— vifi.R.2  ctn20  - 2 ctn  0 — 

E dt2  90 

■ '(A  e>PR  )I5P  0tn  9 - [ 1+v]t-°tn  9 H + V CSC20 


- JlJZ  + (l-  v)R2  9s  _ 0 

902  Eh  30 


(2.13) 


A similar  procedure2  is  used  to  arrive  at  an  equation  for  w: 

*2._  /-<  /i . - . \ _ *->2  _2 

E lit2 


2W  - ±2  R2(s  . p h £h)  - 2<l4v>P^  <L2  + Cl-v^pR4  afs 

•m.  r *v.  2 V-.  *si  O «0.  ~ , 2 


Eh  at 

_ (i~v2)p2r4  a^w 

e2  at ' 


+ ctn  e ^ + £2  . B1  otn  e as 

9^'  90  902  Eh 


E2h  at£ 

ao 


.££•*&  otn  e -4s - + ^ -fo  = 0. 
Eh  902  E 9t290  E 9t2902 


(2.14) 


If  we  introduce  the  new  variables 


w*  = ^5  v*  = ^5 

R r’ 


2 _ t2  _ Et2. 

T 


s*  = 


Rs 

Eh 


then  (13)  and  (14)  reduce  to  the  simpler  non-dimensional  forms 


2.  The  operation  ^(2.9)  + L?(2.10),  where  L.,  = -(ctn0  + -£-) 

o 9 21  60 

and  L2  = 1 + ctn  0 ■$-  + - l4?;v)pR  9 gives  the 

2 50  r E glves  vne 

desired  result. 


3311*11 


M(v*)  = 2(l-v2)  ~S-  + (l-v2)(l-v)  + 2vv»  + v(l-v)JLZ_ 

2 **  4 *2 

9 v /i  .a  3 v . ^J2.  a ±1.  a 3^r* 


and 


- 2 2_-  _ (1-v)  — ^7  + 2v*  ctn2  0 ++(l-v)ctn<: 
— p St  39^ 


39‘ 


9t' 


- 2 ctn  0 221  - (l-v)ctn  0 3.X*.  + [l+v] 

00  aesT2 

- v*  csc20  + ^-~1=  (1-V2)#^ 


ctn  9 Sg 
30 


392  J 


'39 


(2.15) 


L(w*)  = ctn  0 i&±  + ^a*..)  + 2w*  + -d2-W*  + 34w* 

99  3t230  302  3t2302 

,4, 


. (1_v2)0Zwi  + (1.v)  sV  . 2(1+v) 


3t4 


9t2 


32w* 


= ctn  9 tfr  " (1~v2)  + (1~v)s** 


(2.16) 


III.  Acoustic  Wave  Propagation. 

As  in  the  cylindrical  case  [1]  we  have  the  acoustic 
wave  equation 

Acp  - >^2cpTT  = 0 (3.1) 

where  cp  is  the  velocity  potential,  Acp  is  the  Laplacian  in 

spherical  coordinates  and  X2  = -~S|  = -JL.  (c2  £s  the  acoustic 

t2c2  pc2 

speed  of  the  fluid.)  The  pressure  perturbation  p is  again  given 

by 


P = -P*9  (r ,9  ,t) 

T 


> 

.2 


(3.2) 


p* 


PfR  PfE 
t02  P 

where  pf  is  the  fluid  density.  The  applied  stress  s of  (2,2) 
must  of  course  be  the  same  as  the  acoustic  pressure  p at  the 
surface  of  the  sphere.  We  have  in  fact 
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s* 


PfR 

ph 


IV.  The  Interaction  Puflhlm. 

We  now  pose  the  following  problem.  An  incoming  plane 

wave,  cp  Q,  obeying  (3»1)  impinges  on  an  elastic  spherical  shell 
causing  it  to  vibrate  according  to  (2,15)  and  (2.16).  The 
vibrating  body  acts  as  a scatterer  and,  to  a lesser  extent,  as 
a radiator.  The  outgoing  waves  (scattered  and  radiated)  also 
obeying  (3.1),  in  turn  influence  the  nature  of  the  vibrations. 
At  the  surface  of  the  shell,  the  radial  velocity  of  the  fluid, 
tpp(l,0,'n),  must  be  equal  to  the  radial  velocity  of  the  shell, 
w*(0,v).  We  wish  to  determine  the  motion  of  the  sphere  and  the 
pressure  distribution  associated  with  the  incoming  and  outgoing 
waves. 


P =1 


&(z-vA ) 


The  pressure  associated  with  the  incident  wave  is 
taken  to  be  [l] 

r 

Q0e 

0 

so  that  the  initial  velocity  potential  is 


z < t A 

Z > T /X 


90 


Q0X  &(z-t/\)  Q0X 


p*6 


e 


p*& 

0 


For  simplicity  let 

9 
v* 


9 o = P* 

Qo 


q-q*°  , 


w*  = w > 


_ Qo 

p* 

, Qo 
p* 


0|)°  + 1|>  ) - 


v. 


z < t A 
z > v A. 

- x 

p*  A 


(4.1) 


(4,2) 


(4.3) 


n 
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Then  our  boundary  value  problem  is  defined  by  the  following 
set  of  equations: 

L(W)  = pCxlCljQj^ctn  0 - (1-v2)xt^(1>9)t)  + X^1:0*17) 


+ (i-v)xT(i,e,T)3 

(4.4) 

M(V)  = fS(I-v‘ 

-)^(i,9,^) 

(4.5) 

0 

(4.6) 

=xr 

(1,9,t)  [Boundary  Condition] 

(4.7) 

^ o = " 

— e6(WX)  -h  z<tA 

5 b 

0 z > tA  . 

(4.8) 

These  equations  will  be  more  easily  handled  if  Fourier 
transforms  are  first  introduced  to  eliminate  the  time  dependence# 
Denote  the  transform  of  a function  F by  F.  Then  F and  F are 
related  by: 


F(r,0,q-ia)  = 


+ oo 


-oo 


F(r ? e,T)e”aTe”^^T  dT 


(4.9a) 


F(r,0,x)  = ~ eaT 
2n 


where  a is  any  positive  real  number* 

' + oo  _ , 

F(r,0,r,-ia)e1T‘T  dr] 

J -oo 

or,  letting  q-ia  =£ 

1 n-ia+QO  .r 

F(r,0,T)  = FCrjQj^Je1^  dC. 

2lt  J -In -oo 


(4.9b) 


This  generalized  definition  of  the  Fourier  transform  has  been 

r+co  ■ 

used  because  the  usual  definition,  F(r,0,q)  = F(r, 0,T)e",lT^T  dv, 

J -00 


3,  The  operators  L and  M which  appear  here  are  those  given  in 
(2.1?)  and  (2.16). 
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fails  in  the  case  of  the  function  x.^ 

Applying  (4.9a)  to  (3.4)  through  (3.8)  we  obtain 

L(W)  = ip?  C x * ctn  0 + x"(l,0,?)  + (l-v)x  + ?2(l-v2)x]  (4.10) 

where 


= ctn  9 (-JL  - ?2  9 ) + 2 + ?2  ®L. 

l 90  aF  a02  a02 


- (l-v2)?^  - (1-v)?2  + 2(l+v)?2 


M(V)  = ip?  (l-v2)x* 


(4.11) 


where 


M = - -2(1- v2)?2  + (l-v2)(l-v)?If  + 2v  - v (1-v)?2 
2 

- 2 ~ + (l-v)?2  + 2 ctn2  9 - (1-v)?2  ctn2 9 

a0^- 

- 2 ctn  0 + (1-v)?2  ctn  0 -J-  + (l+v)(ctn0  -9 

uu  90 


- esc2  9 + 


90 


90‘ 


A \|)  + X2?2^’  = 0 
i?W(0,?)  = xr(l,0,O 
-e"iXz^ 

^ = i?(&A+i?)  = "f(^)e 


-iXr?  cos  0 


(4.12) 

(4.13) 

(4.14) 


It  can  be  shown,  using  the  method  of  separation  of 
variables  on  (4.12),  that  \jT  is  of  the  form 

r ^ t Jn+f^r^  + CYn+^.A?r)  Pn(cos  9)]. 


4,  It  should  be  noted  that  for  the  functions  we  are  interested 
in, > the  tendency  of  the  integrand  in  (4.9a)  to  become 
infinite  for  large  negative  t is  only  apparent.  Actually, 

- 0 for  -oo  < t /X  < z$  and  W,  V,  and  are  zero  until 
the  wave  hits  the  shell,  i.e.  for  -oo  < VX<-1,  so  that 
integral  of  (4.9a)  always  exists. 
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Since  the  waves  associated  with  \jT  must  be  outwardly  moving,  take 
C = -i  and  write: 

f = ? An(Oh^2)(\Cr)Pn(cos  9) 

where 


(4.15) 


hn2>(^r)  = fe 


(1/2 


[Jwi(X5r)  - 


rJ)Q  may  be  expanded  similarly: 

oo 


uu  n 

i|>°(r,e,0  = -f(0  £ (2nfl)(-i)  X (Hr)P  (cos  0)  (4.16) 

n=0  n 


oo 


= F Bn(OJn(XCr)Pn(cos  0) 
n=0  n n 


where 


/ , r V 1/2 


L(Hr)  = 


From  the  work  of  Lamb  [4]  we  know  that  the  complete  solutions 


for  W and  V may  be  written: 

55(6,0  = 2 W_(C)P  (cos  0) 

m=0  n n 

00  1 


(4.17) 


00 


V(0,o  = - r V (OP  (cos  0)  = -£  W P1  (cos  0 )sin  0 • (4.1® 

' r»  n n n 


n=0 


n ' n 


n=0 


Substituting  these  expressions  into  (4.10),  (4,11)  and 
(4.13),  we  get  three  algebraic  equations  for  the  three  unknowns 
Wn,  Vn  and  An. 


[(l-?2)(n2+n)  - 2 + (l-v2)^  + (l-v)52  - 2(l+v)£2]wn 

= -iKC(l-v)  + (1-v2)£2  - (n2+n)][Anh^2)(\0  + B^XO]  (4.19) 
[2(1+ vK2  - (1-v2)^  + v£2  + 1 + (n+n2-l)(^2-l)]vn 
= -ipOl+v)[Anh^2)(XO  + Bnjn(XO] 


(4.20) 
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\ZBn(K)jn  (W  +xan(0^2)'(XO  = UWn  (4.21) 


These  are  readily  solved  to  give: 
(2nH)(-i)nC? 

Wr  “ xV(i02(r  + iOCcAh(2^'  + iC2h<2;3 
a *■  n ^ n 


(4.22) 


Vn=  - 


An  = “ 


(2n+i)(-i)np(l+v) 


H(U)(£  + iOtc-jXh^  ' + *C2h<2>] 

l[1  + TTlzT' L%T  ]-L(2rwl)  t-i)n 

2 ciXhi2J  + iC0h(2'  K2 
x n 2 n 


(4.23) 


(4.24) 


where 

Cn 


cl  = 


(l-£2)  (n2+n)  - 2 + (1-v2)?4  + (l-v)£2  - 2(l+v)£2 
iP5[(l-v)  + (l-v2)^2  - (n2+n)  ]. 


The  quantities  of  physical  interest  to  us  are  the  stresses  and 
radial  acceleration  associated  with  each  vibrational  mode,  and 
the  total  pressure  distribution.  These  can  all  be  found  at  least 
in  principle  from  (4.22),  (4.23),  (4.24),  (4.15)  and  the  in- 
version formula  (4,9b). 


V.  Numerical  results  for  the  shell. 

The  transforms  of  the  stress  components  will  all  be 
linear  combinations  of  \ and  Vn  (in  general,  gx(0)Wn  + ggO)^), 
and  the  transforms  of  the  radial  accelerations  will  be  -S2Wn. 
From  (4,22)  and  (4,23)  we  see  that  these  expressions  are  regular 
in  ^ except  for  a finite  number  of  poles.  The  theorem  of  resi- 
dues can  therefore  be  used  to  evaluate  the  integral  of  (4,9b). 
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The  computation  will  be  carried  out  in  detail  for  the  three  lowest 
modes.  The  parameters  will  be  taken  as 

p = 17,  6=0,  X2  = 12,  v = .3 


the  values  appropriate  to  a wave  of  infinite  length  impinging  on 
a steel  shell  in  water. 

Zeroth  Mode,  For  this  mode,  the  motion  of  the  sphere 
is  very  simple,  consisting  of  uniform  (9  - independent)  expan- 
sions and  contractions  of  varying  amplitude  and  Deriod, 

The  stresses,  and  Nq,  are  given  by 

Eh  wo 


N = Nq  = 

9 9 1-v  R 


— °R-—  W . 
(l-v)p  0 


w - i -iqXr 

® OlT  4 /V  1 1,2 1 D 4v*  27  4 


2*  -iaX-oo  2C  (C 3“k2C  -I83p+  ik2 ) 

p 2\2 

Here  k2  = ; C = H . 

The  integrand  is  regular  except  for  poles  at 


C = 0,  16.581,  +1  + ,71i 

with  residues  respectively  of 

1,  .0060e-l6*58(X+1J,  and  -1.28l8e-'71(K  +1>sln(S  + 1.89). 

A. 

Jordan's  lemma  can  be  adapted  to  give 

W = 1.00  - 1.28l8e“,71(I  +1)sin(£  + I.89) 

0 2 ^ 

- 0,0060e"l6*^8(X  + 1^.  (5,1) 

A plot  of  -2Nj/RQo  appears  in  Figure  2,  ^ 

5,  These  scales  will  also  be  used  inplotting  the  1st  and  2nd 
modes.  The  reason  for  their  choice  will  become  clear  on 
page  18, 
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The  radial  acceleration  is  found  directly  from  (5.1) 


to  be 


8S,  _ Q0r2  8\  _ Q0R2 

dx2  Ehp  4Eh 


"•71(2.  +1) 

0 K 


+ 0.2697  cos + 1.89) 


0.09^1  sin  + 1*89> 

-16.58(5.  +1)1 
- 0,2427  x k 


4Eh  a^w 
The  plot  of  - * — *■ 

Q0R2  &c2J 


(5.2) 


appears  in  Figure  5* 


There  is  an  alternative  to  the -computational  scheme  we 
have  used,  in  which  the  integrand  is  expanded  in  a power  series 
about  infinity  and  the  series  integrated  term  by  term.  Since 
there  was  a possibility  that  it  might  prove  simpler,  this  was 
tried  out,  the  zeroth  mode  acceleration  being  taken  as  a test 
case. 


02W 


-3a\+  oo 


& iC<?  +1  )Ay 
o l_  P PCe  \ dC 

0T2  2H 


C^-k2  C-l8i  C2+  ik2 


4Ehw 


- ia\-  oo 


TTC 


Q R 


2" 


. . Jl  ( e1^15 


2tc 


Jo 


03) 


C2  + l8iC“3 


- 299 Ck-  49 57i£"^  + 82201C"5 6 *  + ... 


dC 


(f  +1)  - 36(J  +1)2  + |(299)(^  +D3  + ... 


(5.3) 


The  results,  which  appear  in  Table  1,  demonstrate  the  impracti- 
cality  of  this  new  procedure.  It  is  apparent  that  at  least  12 
terms  are  needed  to  find  only  the  first  maximum  to  within  10$ 


5 These  scales  will  also  be  used  in  plotting  the  1st  and  2nd 

modes.  The  reason  for  their  choice  will  become  clear  on 

page  18, 
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First  Mode 


The  stress  components  here  are 

s N = W 1 3v  _ El 

9 9 i -v  n r ao  , n.i 


W . 

1 

3v 

Eh 

T 

R 

R 

"Jl 

(1-vd 

= . AR  ..fw  + M = _ML_  [ w . v" 

(i-v)pL  d-v)p  L 1 \ 


wrvi 


(l-v)pL  aej1  (i-v)p 
3p  -ict\+oo  2e&\  dC 

2ni  -iocX-oo  2i£lf  + ''^ 


J1 

cos  9, 


2iC4  + 3%3  - H3iC2  - 575t  + 575 i 


The  integrand  has  poles  at 


C=  1.17i,  I6.?8i,  and  +3.81  + 0.645i. 

On  computing  the  residues  we  obtain 

wrvx  . 1.9666e-1-17(V  +1)  - O.W*'**  +1) 

-1.8488e-°’625(\  +1)sin  3.811  (J  +1)  + 1.479  . (5.! 

Figure  3 shows  a plot  of  -2N  /RQq  at  the  point  9 = it , 

The  radial  acceleration  is  found  most  easily  by  in- 
verting the  transform,  -^^W^  cos  9. 


? 2 

QQir  a w 
Ehp  "s? 


QoE2  82W1 

“7”  ' ’ o cos  9 

Eh  (3 


-iaX+oo 


71  1 


-iaX-co 


(102C2-1275)e1C(r  +1)dC 

2iCl++38C3-ll3iC2-575C+575i 


= 1.064e“°‘62^\  +1)sin  2.72+3.8l09(^  +1)] 

-3.5335e"l6’58(I  +1)  + 3.H71e‘,1,17(X +1)  (5.5) 
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qIP  9t2J 


is  plotted  for  9 = 7t  in  Figure  6. 


Second  Mode, 

In  general,  the  two  tangential  stress  components  are 
not  equal  to  each  other  for  this  mode.  We  have,  in  fact, 


Ne  = 


Eh 

r 

rT  Eh 

— — 

1-v2 

Le® + v e<p 

j 

N - — 5- 

9 1-v2 

em  + (i)£Q 
cp  y _ 

where 


_ 1 
6 9 5 


, 9V 
w + —r 

- 3£± 

w + 

si 

99_ 

2 Ehp 

_ 

ae_ 

2 

Ehp  _ 

3 cos2  9(-2V2  + iW2) 
+ (3V2  - ^W2) 


_ 1 
£cp  R 


- — 1 

Q0R 

w + v cot  9 

= 

W + V cot  9 

— ■*- 

2 Eh§ 

- mm 

QpR 
EhpL 


3 cos2  0(-V2  + iW2)  - iWr 


1 "iaXt00  [676.8125C3  - 34491.40625Cleit:(l +1)dC 

2 = S J 5 

-iaX-oo 


2 

where 
D 


~iaX+oo  i/ XT') 

= „ A r 8460.1562 5Ce  ^X  ; dC 

Pir  w 

-iaX-oo 


D 


= -.91^^+19.11  iC6+ 123. 725C 1119.475  iC^-3292.5C3+ 

+4212.25  iC2+2756.25C-2756.25i  (5,6) 

The  denominators  have  seven  zeros  5 

16.54  i,  +7*35  + 0.43 i,  +0.86  + 1.79  i,  +0.92  + 0.007  i. 


The  theorem  of  residues  gives 
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w, 


[7.3483(^  + 1)J[+  2 

r _ 3 n 


u 

T 


2e“1*79(t  +1)j12,2253  sin fo, 8578 


2 = O.2130-16*^1'^'  +1)-  2e“0,43(X  + 1)j  0.0426  [sin  7. 3483 (J  +1)J 

- 0,0254  cos 
*• 

+ 1)  H 

0.9242(^  + 1)J  + 2.9634  00s  0.9242(J  + 1) 


2.8318  00s  0.8?79(X  +1)]4  2e",007(\  + 14  3.614? 


sm 


(5.7) 


v2  = 0.0082e"l6*54(X +1)-  2e0,43(x  +1)]  0.0102  sin  7.35(j~ +D 


+ 0.0888  cos 


7.35(£  +1) 


^ 2e-1‘79(X+1)J2. 


sin  0.8578(^ 


|2.8117 

1)]  + 0.8242  cos  0.8579CJ  + 1)]  j-  2e“*007(I  + ^jo,9015 
• sin  jo. 9242 +1)  ] + 0.7395  cos[o.9242(2 +1)] 


(5.8) 


At  0 = it  where  the  greatest  stresses  occur,  NQ  = = 

Q0R  2N 

(Wo-3V0),  A plot  of 2.  f0r  q = it  appears  in  Figure  4, 

(l-v)p  ^ 2 Q0R 

The  radial  acceleration  is  given  by 


d2w 

*2 


_ Q0R2  a2w 
2 ’ Ehp  3t2J 


_ Q0R2  (3  cos  20tl 
Eh  p dv2  l 4 


a2Wc 


fc2  2n\2  -iaX-00 
where  D is  given  in  (5.6). 


1 "*a|+°°  [676,8l25C^-34491.4o625C3]fii^^  +1^C 


D 


£pW. 


0t£ 

+ 


- = 6,6576  e-l6^4(I+1)  + e-°-43(?+1)j  0.564  sin 


7.3483 (~  + 


1)]  - 0.2457  cos [7.3483 (£  +1)  |+  e’1,79^1^  8.939  sin[o.3578 
(£+l)]~  6,9956  cos[o.8578(j£  +1)  |+  e-0.007(^ +1) 

-[0,6982  sin[o.9242(£  +1)]+  0,5966  cos[o,9242(^  +1)]  k (5.9) 
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4Eh 

«or2 


is 


plotted  for  9 = % 


in  Figure  7. 


If  we  examine  Figures  2,  3,  and  4 it  becomes  evident 
that  the  stresses  in  the  lowest  mode  are  very  much  greater  than 

those  of  the  1st  and  2nd  modes,  and  very  likely  those  of  the 

higher  modes  as  well.  In  this  connection  it  should  be  recalled 
that  the  plots  in  Figures  3 and  4,  which  were  made  for  0 = it , 

show  the  largest  stresses  which  can  occur  in  the  first  and  second 

modes  at  any  time.  We  may  therefore  consider  the  resultant 
stresses  in  the  sphere  to  be  predominantly  those  of  the  zeroth 
mode. 


The  case  of  the  acceleration  is  not  so  simple.  A 

picture  of  the  total  acceleration  cannot  be  gained  by  looking 

at  the  lower  modes.  In  fact,  it  would  seem  from  Figures  5?  6, 

oo  2 

and  7,  that  the  series,  £ does  not  represent  the  total 

n=0  dx2 

acceleration  for  all  t . The  rate  of  change  of  radial  accelera- 
tion is  very  great  for  the  first  three  modes  at  t/X  = -1,  9 = it  5 
probably  the  total  radial  acceleration  as  summed  mode  by  mode  will 


00  g 2-y 

be  discontinuous  at  x = 9 = u , This  indicates  that  2 — 

2 n=0  9t* 

does  not  converge  uniformly,  and  therefore  that  ~ ® i^n 


dx‘ 


n=0 


This  difficulty,  associated  with  the  use  of  the  series  expansion 
as  a method  of  solution,  will  be  encountered  again  in  Section  VI 
when  we  discuss  the  resultant  pressure  distribution.  By  antici- 
pating the  results  of  that  section,  we  can  find  the  initial  value 
g^w 

of  -2-i  at  9 = % . 
dx2 


Equation  (2,2)  gives,  for  % A = -1 


n 


-17- 


PhS^  = £Lh  efw 
5tt2  R2p  9t:2 


Qn 


*+Eh  S w 


Q, 


R2  Tt? 


Taking 


we  have 


S = -2Q 


( (6.9) » 9 = 7i ) 


^Eh 


73  =8- 


Qqr2  St 

Thus,  while  we  can  obtain  a very  good  approximation  of  the 
stresses  by  considering  only  the  lowest  mode,  the  same  is 
definitely  not  true  of  the  acceleration.  The  actual  initial 
acceleration  is  about  40  times  as  great  as  the  maximum  acceler- 
ation in  the  zeroth  mode.  The  agreement  obtained  by  considering 
the  first  and  second  modes  along  with  the  zeroth  is  not  appreci- 
ably better;  the  results  still  differ  by  more  than  a factor  of  6. 


Quasi-Static  Case , 

It  is  of  interest  to  compare  the  hoop  stresses  we  have 
obtained,  as  represented  by  those  of  the  zeroth  mode,  with  the 
stresses  for  the  quasi-static  case;  i.e.  for  the  case  in  which 
the  sphere  is  taken  to  be  rigid  and  scattering  is  neglected. 

Since  the  effects  of  only  the  incident  wave  are  con- 
sidered, we  have  for  the  pressure  at  time  t 


Q027tR(tc+R)  Q0(%-  +1) 
~ 4nR2  = 2 


pressure 


= 3 


Q, 


-1  < ^ < l. 

+ 1 < ^ < 00 . 
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The  stresses  are  given  by 

r 

4 


N9  = % = + 


RQ, 


-1  < i + 1 


+ 1 < ~ < oo 


. £?k 


HQ, 


appears  as  the  dotted  line  in  Figure  2,  where  it  can  be 


seen  that  the  very  simple  quasi-static  case  approximates  the 
more  exact  dynamic  case  very  closely.  The  stress  developments, 
while  slightly  out  of  phase,  are  essentially  parallel  with  a 
difference  in  maxima  of  only  1%, 

For  -1  t/X  £ 1,  there  will  be  an  unbalanced  force, 
due  to  the  incident  wave,  acting  on  the  sphere.  This  will  re- 
sult in  an  acceleration  of  the  rigid  body  in  the  z-direction. 

The  magnitude  of  the  force  is 
r % f 2it 


0, 


0 


R2Q0  cos  © sin  0 dcp  = Q mR2  sin20 


= Qqx[  R2  - (tc)2] 

= Q *R2[l  - (t:A)2]. 


Therefore 


lHiR2hp  = q xR2[l  - (r-)23 

9m2  0 X 

4Eh  a2z 


Rc 


*2  * «o[1  - ^ 


4 Eh  9 z tt  a \ 2 i 

~ — 2 = L 1 - J • 

Q0R2  0r  K 

This  has  been  plotted  for  purposes  of  comparison  in  Figures  5>  6, 
and  7 (dotted  curves). 
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The  maximum  acceleration  in  the  zeroth  mode  is  only- 
one  fifth  that  of  the  rigid  sphere.  The  resemblance  is  greater 
for  the  first  mode,  where  the  two  significant  maxima  occur  within 
the  time  | ^|<1,  the  larger  being  one  half  that  of  the  rigid 
sphere.  In  the  second  mode  where  we  begin  to  have  large  nega- 
tive accelerations,  there  is  considerable  difference  in  form 
between  the  dynamic  and  quasi-static  cases  although  the  maxi- 
mum positive  acceleration  of  the  former,  .7,  has  moved  still 
closer  to  the  rigid  body  value  of  1, 

This  value  of  1 is,  we  recall,  very  much  smaller  than 
the  maximum  of  the  total  radial  acceleration  for  the  dynamic 
case  (p.  17). 


VI.  Resultant  pressure  distribution • 

From  (3.2),  (4.2),  (4.15),  (4.24)  and  (4.9b)  the  total 
pressure  is  known  to  bes 


ptotal  = Incident)  + Pn  = -P*<PT  = -Q0X 

-in+  oo 

I 


= Q . 

* o 2 % 


e1^  d£ 


-in  -oo 
-ia+co 


' + ic2h<1) 

C,\h(2)'  + c 

in  2 n 


P (cos  0)d£ 
n 


(6.1) 


where  C-j.  and  C2  are  given  on  p.  10 . 

Unfortunately,  the  order  of  summation  and  integration 
in  (6.1)  cannot  be  interchanged,  i.e.  the  total  pressure  cannot 
be  found  as  the  sum  of  the  pressures  associated  with  the 


Bll-11 


20 


individual  vibrational  inodes.  This  can  be  seen  for  the  specific 
case  r = 1,  x = -X,  and  0 = % as  follows.  At  the  moment  of 

impact,  x = -X,  we  have  for  each  mode  and  for  all  0,  v„  = w 

2*  o * n 

» -u ~ = 0,  Our  equations  of  equilibrium,  (2.9)  and  (2,10) 
St^*  0 14" 

require  that  the  pressure  at  the  surface  of  the  sphere  must  like- 
wise vanish  for  each  mode  at  x = -X,  so  that  the  total  pressure 
would  be  zero  for  all  0,  We  should  expect  however,  from  what 
is  known  of  the  theory  of  scattering  of  plane  waves,  that  the 
pressure  Qq  would  be  doubled  for  9 = % and  not  reduced  to  zero. 
Alternatively,  we  recall  from  p.  16  that  the  total 
radial  acceleration  is  discontinuous  at  x - -X,  9 = it.  Therefore 
the  total  pressure  as  summed  mode  by  mode  will  be  discontinuous, 
indicating  the  non-uniformity  of  convergence  of  the  series  in 
(6,1). 

Because  of  this  peculiarity  in  convergence,  it  is  not 

possible  to  obtain  an  approximate  solution  for  the  pressure  by 

considering  just  the  first  few  terms  of  the  series.  ijh,  (4,15), 

must  be  found  in  closed  form  if  P,  , _ is  to  be  evaluated.  This 

total 

has  not  as  yet  proved  feasible  because  of  the  complexity  of  the 
summation  which  must  be  made.  It  was  noted  however  that  the 
expansions  for  \J>°  and  \jj  are  very  similar  for  r = 1,  and  £ very 

qr 

large  or  (^  - cos  9)  very  small.  This  fact  can  be  used  to  ob- 
tain the  pressure  at  the  surface  of  the  sphere  for  x/\z  cos  9, 

By  comparison  with  the  expression  for  Up  on  page  9 
it  is  seen  that  the  pressure  associated  with  the  incident  wave 
may  be  written 
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PI  = Q = - 


-iaX+oo 

^oi  f iC^/X  oo 

— , J — 206jn(X^r)Pn(cOS  e>(2nH)(-t)n  • 


-iaX-co 


For  C »P  and  r = 1,  Pjj  (cf.  6,1)  becomes,  to  first  order 
<So‘laY  “glC’A® 

PII  = £ _t  £°  -(2wl)(-i)nPn(cos  ©Mt 


s 2 

£ n=0 


and  Pj  reduces  to 


(6.2) 


-iaX+oo 


p Opj  f ® cos(C  - IT  - I 

1 ' " 'z*  -taX-oo  « n=0  C 

Equations  (6,4)  and  (4,l4)  tell  us  that 

*? Y nnn  ^ 


^oi  f piC^A  00  cos(£  - -w-  - S)  n 

~ J -^-p 2 (2rH-l)(-i)nPn(cos  0)d£, 

-daX-oo  ^ n=0  <» 

(6.3) 


C e~^cos  Q » t;  cos(C  - - — ) (2n+l)  (-i)nP  (cos  0) 

C-^OO  n=0  2 2 n 


Therefore,  take 

(K.  - ^o-iK-*/2)003  e » ” sin(C  - - 2)(2»1) 

C->  00  n=0 


2 2' 


»(-i)nPn(cos  9) 


so  that  to  first  orders 


Qn  f 

PTX  - . J-°— 

1 eii:(r 

- cos  0) 

i J 

cos  9 

11  2 % t 

c 

*e  2 

aq 

= V* 

|(oos 

9+1) 

for 

■ 

vl 

o 

cos  9) < < ,05 

= 0 

for 

<?  - 

cos  9)  < 0 

= 0 for  (^  - cos  9)  < 0 (6,4) 

Higher  order  terms  may  be  obtained  in  the  same  way. 
To  second  orders 


P = Jjf  “ r(Sln«  - f - f)Cl  - ^ 

11  2nJ  £ n=0  L 2 2 e 


2X~  J n 
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Qo 

2% 


f lC-cA  2, 


e 


slnK  - g*  - |)  f cos(C  - jp  - f)n(nH) 


P - 

I " " 2* 


^ n=0  L 

' 2i2 

fcti  ei(C  -f'V 
z 2 _ J 

C2w-i)(-i)nPnCoos  9)ae  (6.5) 

" I^t/X  oo 

cos(C 

- f5  - |)  since  - ™ - f)n(nH) 

J ^ n=0 

£ 2t:2 

• (2nH)(-i)Pn(cos  9)d£  . 


(6.6) 


In  this  case 
-iCcos  9 


(2n+l)(-i)nP  (cos  9) 


ooso:  - a -|) 


since  - - ^)n(iM-l) 

_2 


therefore  by  analogy,  to  second  order: 


-ICC  - f)oos  9 
e * 


oo  n 

->  Z (2nH)(-i)  P 


sinCC  - f - *) 


C ->  oo  n=0 

cos(£  - - ^)n(rn-l) 

+ 2C2 


[life] 


This  does  not  correspond  exactly  to  the  first  two  terms  in  Pjj, 
therefore  we  must  subtract 
% Q-i(C  - ^)cos  9 

1_ — » ? sin(e  - ?!  - 5)(2"-lH-i)n  P„(cos  9). 

t 2 2 


We  also  have,  again  to  second  order 

<P+1) 


Q-iC  cos  9 ie-i(£  - |^cos  ®i 

'V  T " T'"r  J,  " ™ 


CO 

— *(p+l)  T. 

JZ->oo  n=0 


oos(C--f!-  |) 


i sinCC  - f*  - f) 


(2nH)(-i)nPn(cos  Q). 
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Pjj  can  now  be  written; 

IC'tA 


p _ 

II  " 2tc 


e 


e 


II  2tc 


+ $ + l)(* 
iC(^  ••  cos  9) 


. . -i(C  “ §)cos  0 

-i(C-it/2)cos  9 f 2 

“ 2 C „ 

-i£  cos  9 ie“I(C  - ^)cos  9 


) 


dC 


i § COS  9 p+1 


e 


Tt 


+ — (1  + ie  2 

% 

2C 


i --  cos  9 


) 


1C  ei|  cos  9 


A i §(cos  9+1)  /T 

= Q0  [e  2 - (^ 

Tt  /c  . , i | cos  9 

+ - cos  9)e  2 

= 0 


cos  9)  (1+  p)  (1  + e1  2(cos  9 +1) 

for  0 < (r  - cos  9)  « .05 

for  - cos  9)  < 0,  (6.7) 


Additional  terms  will  be  of  little  value  since  our 

expansion  is  valid  only  for  £ »(3,  ox1  (t/\  - cos  9)  < < .05. 

The  terms  we  have  found  so  far  however  are  sufficient  to  tell 

us  some  things  of  importance. 

The  incoming  wave  will  reach  the  point  (1,9)  on  the 

sphere  at  time  t/X  = cos  9,  The  initial  pressure  for  each  9 

6 

is  given  by 


6.  The  elastic  waves  in  the  shell  will  travel  more  rapidly -than 
the  acoustic  waves  and  will  result  in  a pressure,  P 4 0,  at 
(1,9)  before  the  time  t/X=  cos  9,  9 However,  this 

effect  is  negligibly  small  compared  with  the  one  we  are 
considering. 
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Ptotal  = pi  + RePn 

= Q + Q jcos(?Kcos  9+1))  - l8[f  - cos  9][l+cos(~(cos9+l))] 

O O l C.  N 1 C 

+ ~[^  - COS  Q][cos(|  cos  0)]^  . (6.8) 

= + Q cos  [ -(cos  9+1)]  • (6.9) 

o o 2 

At  9 = it,  the  outermost  point  of  the  sphere,  the  first 
effect  is  that  of  a plane  wave  hitting  a rigid  wall  so  that  we 
have 

PII  = plS  Ptotal  = 2V 

At  9 =ic/2,  the  wave  just  grazes  the  sphere  and  therefore 

PII  = °5  Ptotal  = V 

As  9 varies  from  tt  to  it/2,  the  initial  pressure  varies  contin- 
uously from  2Qq  to  Qq. 

If  the  sphere  were  rigid,  the  steady  state  pressure 
distribution  \>rould  be  given  by 

Ptotal  = °<  e<  * . [?] 

The  results  of  section  V indicate  that  we  will  have  asymptotic 
values  of  RQ0/2  for  the  stresses  and  zero  for  the  radial  ac- 
celeration, which  also  correspond  to  a uniform  pressure  of 


At  present,  this  is  about  all  that  can  be  said  on  the 
subject  of  the  pressure  distribution.  A completely  satisfactory 
way  of  dealing  with  the  problem  will  not  be  had  until  it  is 
possible  to  find  the  sum  in  (6.1). 
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VII.  Conclusions. 

Both  the  stress  and  the  acceleration  for  each  mode  can 
readily  be  computed  to  any  desired  accuracy.  For  the  stress,  the 
zeroth  mode  is  by  for  the  most  important  and  this  is  closely 
approximated  by  the  zeroth  mode  of  a siinple  quasi-static  system 
(p.  18).  The  acceleration  on  impact  of  the  outermost  portion 
of  the  shell  (9  = it)  can  be  found  exactly  and  is  seen  to  differ 
markedly  from  the  accelerations  associated  with  the  individual 
vibrational  modes. 

The  problem  here  considered,  apart  from  its  intrinsic 
interest,  should  serve  as  a valuable  guide  in  the  solution  of 
similar  problems  involving  obstacles  of  more  complicated  shape. 
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.pproximating  the  zeroth  mode  acceleration  by  a power  series  expansion, 


1.125 


SI 


Hoop  stress  for  the  second  mode  of  9 


■H-M  32 


Fig.  5.  Radial  deceleration  for  the  zeroth  mode. 


Fig. 7.  Radiol  acceleration  for  the  second  mode  at  9*r 
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Branch  Office 
346  Broadway 

New  York  13,  New  York  (1) 


Commanding  Officer 

Office  of  Naval  Research 

Branch  Office 

1000  Geary  Street 

San  Francisco,  California  (1) 

Commanding  Officer 

Office  of  Naval  Research 

Branch  Office 

1030  Green  Street 

Pasadena,  California  (1) 

Officer  in  Charge 
Office  of  Naval  Research 
Branch  Office,  London 
Navy  No.  100 

FPO,  New  York,  N.  Y.  (5) 

Library  of  Congress 
Washington  25,  D.  C. 

Attn:  Navy  Research  Section  (2) 

Commanding  Officer 

Office  of  Naval  Research 

Branch  Office 

844  N.  Rush  Street 

Chicago  11,  Illinois  (1) 


II:  Department  of  Defense  and  other  Interested  Government  Activities 


a)  General 

Research  & Development  Board 
Department  of  Defense 
Pentagon  Building 
Washington  25*  D.  C. 

Attn:  Library(Code  3D- 1075)  (1) 

Armed  Forces  Special  Weapons 
Project 

P.  0.  Box  2610 
Washington,  D.  C. 

Attn:  Lt.  Col.  G.  F.  Blunda  (2) 

Joint  Task  Force  3 

12  St.  & Const.  Ave,  , N.W# 

(Temp.  U) 

Washington  25,  D.  C. 

Attn:  Major  B,  D.  Jones  (1) 


b)  Army 

Chief  of  Staff 
Department  of  the  Army 
Research  & Development  Division 
Washington  25,  D.  C. 

Attn:  Chief  of  Res.  & Dev,  (1) 

Office  of  the  Chief  of  Engineers 
Assistant  Chief  for  Works 
Department  of  the  Army 
Bldg.  T-7,  Gravelly  Point 
Washington  25,  D.  C, 

Attn:  Structural  Branch 

(R.  L.  Bloor)  (1) 

Engineering  Research  & Development 
Laboratory 

Fort  Belvoir,  Virginia 

Attn:  Structures  Branch  (1) 
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Office  of  the  Chief  of  Engineers 
Asst.  Chief  for  Military 
Construe tion 
Department  of  the  Army 
Bldg.  T-3,  Gravelly  Point 
Washington  25,  D.  C, 

Attn:  Structures  Branch 

(M.  F.  Carey)  (1) 

Protective  Construction 
Branch  (I.  0.  ThornleyXl) 

Office  of  the  Chief  of  Engineers 
Asst,  Chief  for  Military 
Operations 

Department  of  the  Army 
Bldg.  T-7,  Gravelly  Point 
Washington  25,  D.  C. 

Attn:  Structures  Development 

Branch  (W.  F.  Woollard)(l) 

U.  S.  Army  Waterways  Experiment 
Station 
P.  0.  Box  631 
Halls  Ferry  Road 
Vicksburg,  Mississippi 
Attn:  Col.  H.  J.  Skidmore  (1) 

The  Commanding  General 
Sandia  Base,  P.  0.  Box  5100 
Albuquerque,  New  Mexico 
Attn:  Col.  Canterbury  (1) 

Operations  Research  Officer 
Department  of  the  Army 
Ft,  Lesley  J.  McNair 
Washington  25,  D.  C. 

Attn:  Howard  Brackney  (1) 

Office  of  Chief  of  Ordnance 
Office  of  Ordnance  Research 
Department  of  the  Army 
The  Pentagon  Annex  #2 
Washington  25,  D.  C. 

Attn:  ORDTB-PS  (1) 

Ballistics  Research  Laboratory 
Aberdeen  Proving  Ground 
Aberdeen,  Maryland 
Attn:  Dr.  C,  W,  Lampson  (1) 

c ) Navy 

Chief  of  Naval  Operations 
Department  of  the  Navy 
Washington  25,  D.  C, 

Attn:  OP-il  (1) 

OP-363  (1) 


Chief  of  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Director  of  Research  (2) 

Code  1+23  (1) 

Code  *f42  (1) 

Code  lf21  (1) 


Director,  David  Taylor  Model  Basin 
Department  of  the  Navy 
Washington  7,  D.  C. 

Attn:  Code  720,  Structures 


Division 

(1) 

Code  7*+0,  Hi- Speed 

Dynamics  Div. 

(1) 

Commanding  Officer 

Underwater  Explosions  Research  Div. 

Code  290 

Norfolk  Naval  Shipyard 
Portsmouth,  Virginia  (1) 

Commander 

Portsmouth  Naval  Shipyard 
Portsmouth,  N.  H, 

Attn:  Design  Division  (l) 

Director,  Materials  Laboratory 
New  York  Naval  Shipyard 
Brooklyn  1,  New  York  (1) 

Chief  of  Bureau  of  Ordnance 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Ad-3,  Technical  Library  (1) 


Rec,  P.  H.  Girouard  (1) 

Naval  Ordnance  Laboratory 
White  Oak,  Maryland 
RFD  1,  Silver  Spring,  Maryland 
Attn:  Mechanics  Division  (1) 

Explosive  Division  (1) 

Mech.  Evaluation  Div.  (1) 


Commander 

U.  S.  Naval  Ordnance  Test  Station 
Inyoke rn.  California 
Post  Office  - China  Lake,  Calif, 
Attn:  Scientific  Officer  (1) 

Naval  Ordnance  Test  Station 
Underwater  Ordnance  Division 
Pasadena,  California 
Attn:  Structures  Division  (1) 
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Chief  of  Bureau  of  Aeronautics 
Department  of  the  Navy 
Washington  25?  D.  C. 

Attn:  TD-41,  Technical  Library 


Chief  of  Bureau  of  Ships 
Department  of  the  Navy 
Washington  25,  D.  C. 

Attn:  Code  P-314  (1) 

Code  C-313  (1) 

Officer  in  Charge 
Naval  Civil  Engr.  Research  & 
Evaluation  Laboratory 
Naval  Station 

Port  Hueneme,  California  (1) 

Superintendent 

U.  S.  Naval  Post  Graduate  School 
Annapolis,  Maryland  (1) 
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d)  Air  Forces 

Commanding  General 
U.  S.  Air  Forces 
The  Pentagon 
Washington  25,  D.  C. 

Attn:  Res.  & Development  Div.  (1) 

Deputy  Chief  of  Staff,  Operations 
Air  Targets  Division 
Headquarters,  U.  S.  Air  Forces 
Washington  25,  D.  C. 

Attn:  AFOIN-T/PV  (1) 

Office  of  Air  Research 
Wright-Patterson  Air  Force  Base 
Dayton,  Ohio 

Attn:  Chief,  Applied  Mechanics 

Group  (1) 

e ) Other  Government  Agencies 

U.  S.  Atomic  Energy  Commission 
Division  of  Research 
Washington,  D.  C.  (1) 

Director,  National  Bureau  of 
Standards 
Washington,  D.  C. 

Attn:  Dr.  W.  H.  Ramberg  (l) 
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Brown  University 
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Stanford  University 

Stanford,  California  1 

Professor  R.  M.  Hermes 

College  of  Engineering 
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Department  of  Aeronautical 

Engineering  & Applied  Mechanics 
Polytechnic  Institute  of  Brooklyn 
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Brown  University 
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Dr«  So  Raynor 

Armour  Research  Foundation 
Illinois  Institute  of  Technology 
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Professor  E#  Reissner 
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Professor  M.  A.  Sadowsky 
Illihois  Institute  of  Technology 
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Professor  V.  L Salerno 

Department  of  Aeronautical  Engineering 

Renssalaer  Polytechnic  Institute 

Troy,  New  York  1 

Professor  M.  G.  Salvador! 

Department  of  Civil  Engineering 

Columbia  University 
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New  York  27,  New  York  1 

Professor  J.  E Stallmeyer 
Talbot  Laboratory 
Department  of  Civil  Engineering 
University  of  Illinois 
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Technology  Center 
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Professor  R.  G.  Sturm 
Purdue  University 

Lafayette,  Indiana  ± 

Professor  F.  K.  Teichmann 

Department  of  Aeronautical  Engineering 
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